Astronomy & Astrophysics manuscript no. santos'print 


© ESO 2010 


November 10, 2010 





The HARPS search for southern extrasolar planets* 
XXV. Results from the metal-poor sample 

q". N.C. Santos 1 ' 2 ' 3 , M. Mayor 3 , X. Bonfils 3 ' 4 X. Dumusque 1 ' 3 , F. Bouchy 5 , P. Figueira 1 , C. Lovis 3 , C. Melo 6 , 
—i ' F. Pepe 3 , D. Queloz 3 , D. Segransan 3 , S. G. Sousa 1 , and S. Udry 3 



, 1 Centro de Astrofisica, Universidade do Porto, Rua das Estrelas, 4150-762 Porto, Portugal 

■ 2 Departamento de Ffsica e Astronomia, Faculdade de Ciencias, Universidade do Porto, Portugal 

, O ' 3 Observatoire de Geneve, Universite de Geneve, 51 ch. des Maillettes, 1290 Sauverny, Switzerland 

. 4 Laboratoire d'Astrophysique, Observatoire de Grenoble, Universite J. Fourier, CNRS (UMR5571), BP 53, F-38041 Grenoble, 

' Cedex 9, France 

5 Institut d'Astrophysique de Paris, UMR7095 CNRS, Universite Pierre & Marie Curie, 98bis Bd Arago, 75014 Paris, France 

i | 6 European Southern Observatory, Casilla 19001, Santiago 19, Chile 

■ Accepted for publications 



ABSTRACT 



o 



> 

o 

(N 



Searching for extrasolar planets around stars of different metallicity may provide strong constraints to the models of planet 
formation and evolution. In this paper we present the overall results of a HARPS (a high-precision spectrograph mostly dedicated 
to deriving precise radial velocities) program to search for planets orbiting a sample of 104 metal-poor stars (selected [Fe/H] 
below —0.5). Radial velocity time series of each star are presented and searched for signals using several statistical diagnostics. 
Stars with detected signals are presented, including 3 attributed to the presence of previously announced giant planets orbiting 
the stars HD171028, HD181720, and HD190984. Several binary stars and at least one case of a coherent signal caused by activity- 
related phenomena are presented. One very promising new, possible giant planet orbiting the star HD 107094 is discussed, and 
the results are analyzed in light of the metallicity-giant planet correlation. We conclude that the frequency of giant planets 
orbiting metal-poor stars may be higher than previously thought, probably reflecting the higher precision of the HARPS survey. 
In the metallicity domain of our sample, we also find evidence that the frequency of planets is a steeply rising function of the 
stellar metal content, as found for higher metallicity stars. 



Key words, planets and satellites: formation - Stars: abundances - planetary systems - Techniques: spectroscopic - Techniques: 
radial velocities 
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1. Introduction 

The past 15 years have seen a plethora of exoplanet discov- 
eries, with more than 450 announcements so faiQ. Mostly 
thanks to the development of the radial-velocity tech- 
nique, some of the discovered planets have been found 
to have masses as l ow as ~2 times the mass of our Earth 
( Mavor et al. 20091 ). The holy grail of exoplanet searches, 



namely the first detection of an exo-Earth, is closer than 
ever. 

The increasing number of detected planets is also pro- 
viding strong cons traints for the models of planet forma- 
tion and evolution ( Udrv fc Santosl[2007 : lldafc Linll2004a[ 



* Based on observations collected at the La Silla Parana 
Observatory, ESO (Chile) with the HARPS spectrograph at 
the 3.6-m telescope (ESO runs ID 72.C-0488, 082.C-0212, and 
085.C-0063). 



See updated table at http://www.exoplanet.eu 



Mordasini et al.ll2009al ). Two main models of giant planet 



formation have been discussed in the literature. On the 
one hand, the core-accretion model tells us that giant 
planets can be formed by the accretion of solids in a pro- 
toplanetary disk, building up a ~10 M m core followed by 



rapid a gglomeration of gas (IPqllack et alj|l996tllda fc Lin 



l2004at iMordasini et alj l2009aj) . On the other hand, the 



"disk instability" model suggests that giant planets may 
be the outcome of direct gravita tional instability of the gas 
(|Bosslll997t iMaver et alTl2002h . The two models predict, 
however, very different outcomes for the planet formation 
process as a function of the chemical com position of the 



disk (i.e. density of solids) and stellar mass ([Matsuo et al 



20071 ). In particular, while the core accretion model pre- 



dicts a strong metallicity-planet correlation and a higher 
preval ence of giant planets around higher mass sola r-type 
stars ([Mordasini et al.ll2009bt lLaughlin et al.l 120041 ) . such 
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trends are l ess evident for the disk instability process (Boss 
20021 l200fih . 

Concerning stellar metallicity, the first large, uniform 
spectroscopic studies comparing large samples of stars 
with an d without plan ets have confirmed former sus- 
picions ( Gonzalez! 1997 ) that stars hosting giant plan- 
ets are more metal-rich (on average) than single , field 



stars (|Santos et al.ll200lL 120045 iFischer fc ValentilliooBI ). 
Stellar mass is also suspected to have a strong influ- 
ence on the planet forma tion process ((Johnson et al 



2007 



Lovis fe Mavorl 120071 ). but a clear understanding 
of this picture is required. This need has inspired the 
construction of specific samples to search for planets 
around stars of different mass /evolution a ry states (e.g 



Bonfil s et al. 2005 ; Setiawan et al. 2005t Johnson et al 



Da Silva et al 
20091) . 



2007t lLovis fc Mavorl 120071) and di ff erent chemica l com 



positions (e.g. Tinney et al. 20031 : Fischer et ahl 20051: 

;t all " 



2006; Santos et all 2007: Sozzctti et al 



The HARPS guaranteed time observations (GTO) pro- 
gram started to follow sever al different sample s of solar- 
type stars in October 2003 (jMavor et al.ll2003l ). The un- 
paralleled long-term precision of HARPS allowed discov- 
ery of several planets among the targets, including the 
large majority of the known planet s with masses near 
the mass of Neptune or below (e.g . ISantos et al.ll2004a ; 



Lovis et alJbood: iMavor et alJliob" 



One of the HARPS GTO sub-samples was built to ex- 
plore how frequently giant planets orbit metal-poor stars 
to determine the metallicity limit below which no giant 
planets can be observed. In this paper we present the over- 
all results of this program. In Sects. [2] and [3] we present 
our initial sample and observations. In Sect.0] we show 
how the sample was cleaned a posteriori of the presence 
of stellar binaries and fast-rotating, active stars. The or- 
bital parameters for some of the binaries in the sample 
are also listed. The remaining stars are then discussed 
in Sect. [5] The three planets already announced and one 
more promising candidate are presented. We conclude in 
Sect.El 

2. The HARPS metal-poor sample 

To explore the low-metallicity tail of the planet-host 
stars' distribution, one of the samples studied within 
the HARPS GTO program followed 104 metal-poor or 
mild metal-poor solar-type stars. This sample was cho- 
se n based on a pre l imina ry version of the large catalog 
of Nordstrom et ahl ([2004). From this catalog, we took all 
late-F, G, and K stars (b - y >0.330) south of +10° of 
declination that have a visual V magnitude brighter than 
12. From these, we then excluded all identified visual and 
spectroscopic binaries, all stars suspected to be giants, 
and all those with measured projected rotational veloc- 
ity v sin? above ~6.0km.s~ 1 (to indirectly exclude the 
most active stars, for which it is not possible to achieve 
a high radial- velocity precision - ISaar fc Donahue 1997 ; 
Santos et al. 2000i Paulson et al. 20021) . Finally, we con- 



sidered only those targets with estimated photometric 
[Fe/H] between —0.5 and —1.5. As seen below, a detailed 
spectroscopic analysis has later shown that some of the 
targets fall significantly out of this metallicity interval (see 
Fig.H). 

The final 104 stars in the sample (Table[T]) have their 
V magnitudes between 5.9 and 10.9, distributed around 
an average value of 8.7. After a 15-minute exposure, these 
magnitudes allow us to obtain an S/N high enough to de- 
rive radial- velocities with a precision better than lms^ 1 
for the majority of the targets. We note, however, that not 
all of these were interesting targets for a high-precision ra- 
dial velocity planet search program (Sect. 2]). 

Precise spectroscopic atmospheric parameters and 
masses for the tar gets were derived in a separate paper 
(jSousa et ahl 120101 ) . In Fig. Q] we plot the distributions of 
metallicities, temperatures, stellar masses, and V magni- 
tudes for the sample. In the upper-left panel, the dashed 
line represents the whole sample, while the filled and dot- 
ted lines represent the distributions for stars with high 
amplitude long term radial velocity trends (discussed in 
Sect. 15.3.12]) and stars with known planets (Sect. l5T3~Tj) . 
respectively. Overall, the metallicity distribution peaks 
around — 0.65dex and has a strong fall above — 0.4dex. 
The sample includes, however, several stars with spec- 
troscopic metallicity above ~— 0.5dex, even though the 
initial cut, based on photometric metallicity estimates, 
was set at — 0.5dex. Some of these outliers can be ex- 
plained by the error bars, both in the spectroscopic and 




20 40 5 10 15 

Number of measurements Weighted RMS [m/s] 

Fig. 2. Left: Histogram of the number of radial- velocity 
measurements for the stars in our sample. Right: 
Distribution of weighted rms values for our stars with at 
least 5 radial-velocity measurements. For clarity, the dis- 
tribution is only shown for rms values below 15 ms -1 . 

Table 1. List of targets, coordinates (equinox 2000.0), 
and magnitudes. Full table available online. 



Star 


alpha 


delta 


mv 


HD224817 


00:00:58.2 


-11:49:25 


8.4 


HD967 


00:14:04.4 


-11:18:41 


8.4 


HD11397 


01:51:40.5 


-16:19:03 


9.0 


HD 16784 


02:40:38.7 


-30:08:07 


8.0 


HD 17548 


02:48:51.8 


-01:30:34 


8.2 
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Fig. 1. Distributions of metallicity (top left), effective temperature (top right), stellar mass (bottom left), and visual 
magnitude (bottom right) for all the stars in our sample. 



(initial) photometric metallicities (the sample had been 
cut using the latter). The most extreme cases (HD 62849, 
HD 123517, HD 144589, and CD-231087, all with [Fe/H] 
above —0.3 dex) are less easy to explain, though. HD 62849 
([Fe/H]=-0.17) is a close visual binary (see Sect. Oj) . We 
cannot exclude light contamination from the secondary 
star having biased the derivation of the stellar parameters. 
In any case, this star was later excluded from the sample 
(Sect.HJ. CD-2310879 is the hottest star in our sample, 
and it presents a high activity level fSect. 15.33)) . We can- 
not reject that the presence of active regions has induced 
a systematic error in our parameters. No clear explana- 
tion exists for the remaining two stars, namely HD 123517 
and HD 144589. Interestingly, both present clear long- 
term trends in the radial-velocity data (see discussion in 



Sect. 15.3.12)) . Their high metallicities may be a signature of 
long period giant planets. 

The mass and effective temperature distributions are 
rather peaked, with maxima around 0.8 M Q and 5700 K. 
Both distributions are roughly Gaussian shaped, although 
the former shows a tail towards higher masses. All stars 
except HD 221580 (Sect 14.21) have spectroscopic parame- 
ters (temperature and surface gr avity) typical of d warfs 



or subgiants (for more details see Sousa et al. 2010l ) 



In Tabled we present the basic stellar parameters for 
the sample stars, excluding objects discussed in Sect.HJ 
The effective temperatures, surface gravities, and stel- 
lar metallicities were taken from the d e tailed spectro- 
scopic analysis presented in Sousa et al. ( 2010l ). The B- 



V values were taken from the Hipparcos catalog (ESA 
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Table 2. Stars excluded from the sample. 



Star 




Comment 




BD- 


00 4234 


Active/Fast rotator 


BD- 


03 2525 


SB2 




CD- 


45 2997 


SBlt 




CD- 


43 6810 


SB1 




HD16784 


SB1 




HD25704 


Close visual binary 


HD62849 


Close visual binary 


HD128575 


Active/Fast rotator 


HD134113 


SBlf 




HD161265 


SB2 




HD164500 


SB2 




HD175179 


Close visual binary 


HD187151 


SB2 




HD197890 


Active/Fast rotator 


HD221580 


Giant 




HD224347 


SBlf 




f Orbital fit was possible. 
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Fig. 3. Radial-velocity measurements of HD16784 and 
CD-436810, two SB1 binaries from the sample. 



19971) . when available, and for the remaining stars from 



the NOMAD catalog (jZacharias et all 120041 ). The esti- 
mate of the projected rotational velocity v sini was de- 
rived using a calibration based on the width of the HARPS 
cross-correlati on function, s i milar t o the one presented for 
CORALIE in ISantos etail (|2002h . Due to different ob- 
serving uncertainties, we adopted a value of <lkms _1 
for stars presenting an estimated v sini below this value. 
The activity level of the star, denoted here by the Smw 
and logi?^^- valuefH, were derived from the analysis of 
the Ca II H and K lines in the HARPS spe ctra, using a 



metho dology similar to the one presented in ISantos et al 
(|2000l) . 



3. Observations 

From October 2003 to July 2010, a total of 1301 precise 
radial- velocity measurements were obtained with HARPS 
for the stars in our program. The observations were done 
as part of the HARPS GTO program, and later comple- 



Both are presented since the latter depends on B-V, a pa- 
rameter that we do not control. 
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Fig. 4. Top: Radial- velocity measurements of HD 134113 
as a function of time, and the best Keplerian fit to the 
data. The residuals of the fit are shown in the lower panel. 
Center and bottom: Same for HD224347 and CD-452997, 
respectively. 



mented by two separate follow-up programs for a few stars 
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Table 3. List of targets with number of measurements, time span between first and last measurements, average radial 
velocity, stellar atmospheric parameters, B-V color, activity level, and v sin?. Full table available online. 



Star 


N 


Span 


<RV> 




log g 


[Fe/H] 


Smw 


c(Smw) 


B-V 


log R' HK 


v sini 






[days] 


[krns" 1 ] 


[K] 


[cgs] 












[kms -1 ] 


HD967 


34 


2132 


-24.3 


5568±17 


4.53±0.02 


-0.68±0.01 


0.176 


0.003 


0.645 


-4.91 


<1 


HD11397 


33 


2075 


41.1 


5564±26 


4.46±0.04 


-0.54±0.02 


0.184 


0.002 


0.693 


-4.89 


<1 


HD17548 


10 


1504 


-15.4 


6011±26 


4.44±0.02 


-0.53±0.02 


0.162 


0.005 


0.529 


-4.91 


<1 


HD17865 


21 


1853 


31.9 


5877±24 


4.32±0.03 


-0.57±0.02 


0.159 


0.002 


0.568 


-4.96 


<1 


HD22879 


36 


1857 


120.4 


5884±33 


4.52±0.03 


-0.81±0.02 


0.165 


0.003 


0.554 


-4.91 


<1 



(ESO programs ID 082.C-0212 and 085.C-006^|) aimed 
at searching for Neptune-like planets orbiting a sample of 
moderately metal-poor stars. Most of the stars have 5 or 
more measurements (Fig.[2]). 

During the first three years of the survey, the mea- 
surements were obtained with ~2ms _1 precision. The de- 
fined exposure times were not long enough to average out 



the " noise" coming from stellar oscillations (|Bouchv et al 



20051) . Since 2006, this strategy has been revised, and the 
full precision of HARPS was exploited. Exposure times 
were then set to be 900 s, allowing the precision of the in- 
dividual measurements to be increased to below lms . 
The observing strategy used does not avoid other sources 
of noise, includi ng those coming from stellar granulation 
and activity (e.g. lDumusaue et al.ll2010l) . As seen in Fig.[21 
if we consider only objects with at least 5 measurements, 
most of the stars present a radial- velocity rms of ~1- 
2.5ms" 1 , as expected. In this figure we limited the rms 
range to values up to ~15 ms" 1 for the sake of clarity. 
Only 10 objects fall above this limit. 

The data were reduced using the latest version of the 
HARPS pipeline. Thi s includes a correction fo r the secu- 
lar accelleration (e.g. Zechmeister et al. 20091) . for which 
we u sed Hipparcos parallaxes and proper motions (ESA 
1997), whenever available. 



4. Cleaning the sample 

A few stars from our original sample were found to be 
unsuitable targets for a planet search survey for different 
reasons (binarity, activity, high rotation - Tabled . In this 
section we present these targets, which will be excluded 
from the rest of the discussion. 



4.1. Binaries 

Though known binary stars were a priori excluded from 
the sample, a few SB1, SB2, or close visual binary systems 
passed through our first selection criteria. After a few mea- 
surements, these cases were usually discarded from the 
sample. 

Three of these stars, HD25704 (number of measure- 
ments, Nmes=20), HD62849 (Nmes=17), and HD175179 



In this latter case only for a star discussed in 15.3.41 



(Nmes=3), were found to be close visual binary sys- 
tems, after visual inspection of the guiding camera images. 
These have been at some point excluded from the sample 
given that the light contamination from the stellar com- 
panion in the HARPS fiber precludes any determination of 
precise radial velocities. This contamination is seeing de- 
pendent. These stars are not adequate for a planet search 
survey. 

BD-03 2525 (Nmes=2), HD161265 (Nmes=2), 
HD164500 (Nmes=2), and HD187151 (Nmes=l) were 
also excluded from the sample since an inspection of the 
cross-correlation function revealed that they are SB2 bi- 
naries. B P— 03 2525 was a l so cla ssified as a spectroscopic 
binary in IPourbaix et al. (|2004l) (after the be ginning of 
our survey). 

Five other stars in our sample were found to be SB1 
binaries. For three of these (see next section), we ob- 
tained enough HARPS points to obtain a reliable orbital 
solution fSect. l4Tl~l"j) . For the remaining two, HD16784 
(Nmes=3; span=272 days) and CD-43 6810 (Nmes=9; 
span=1157days), the reduced number of measurements or 
their span did not allow us to fit any Keplerian function 
to the data (Fig.©. 

4.1.1. Binaries with orbital fits 

For three of the binary stars in our sample, we gathered 
enough measurements to allow derivation of at least a ten- 
tative orbital solution. 

HD 134113: Though previously known to be a spectro- 
scopic binary ([Latham et al] 120021 ). we missed this clas- 
sification when preparing the sample. A total of 28 pre- 
cise radial-velocities were thus obtained for this star. A 
Keplerian fit to the data confirms and refines the previ- 
ously suspected solution, with P= 201.7days, e=0.89, and 
K=3.9kms~ 1 . The complete set of orbital parameters is 
listed in TablelH Consider ing a stellar mass of 0.88 solar 
masses ( Sousa et al. 2010h . this signal is compatible with 
a minimum mass companion with 48 times the mass of 
Jupiter, a possible brown dwarf companion. 

HD 22434 7: There is no previous reference in the litera- 
ture of this star being a spectroscopic binary. Our 8 precise 
radial velocities obtained with HARPS show, however, a 
clear high- amplitude signal. The radial-velocities are well 
fit with a Keplerian function with period 6380 days, eccen- 
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Table 4. Elements of the fitted orbits for SB1 binaries. 



HP 134113 HP 224347 CP-45 2997t 



p 


201.674±0.008 


6380±2500 


22.6 


[d] 


T 


2 453 984.87±0.01 


2 453 665±78 


2 453176.8 


[d] 


a 


0.66 


7.0 


0.15 


[AU] 


e 


0.888±0.001 


0.56±0.11 


0.50 




V r 


-60.19±0.06 


1.38±0.17 


39.2 


[kms" 1 ] 


UJ 


163.1±0.2 


221±4 


172 


[degr] 


K\ 


3907±42 


2890±210 


9500 


[ms- 1 ] 


AM 


121,138 10~ 6 


9 062.746 10" 6 


1301.588 10" 6 


[Mq] 


a{0 - C) 


1.55 


0.56 




[ms- 1 ] 


N 


28 


8 


6 




m,2 sini 


48 


210 


92 


[M Jup ] 



f Number of points is similar to number of free parameters; no errors available. 



tricity of 0.56, and K=2.9kms _1 . The large uncertain- 
ties (see Tabled come from our data not covering one 
full orbital period, and the phase space is not covered 
well. This orbital solution should be considered as pre- 
liminary. Consider ing a stellar mass of 0.95 solar masses 
(jSousa et al.ll2010h , this signal is compatible with the pres- 
ence of a minimum mass companion with 210 times the 
mass of Jupiter orbiting HD 224347. 

CD— 45 2997: 6 radial- velocity measurements of this 
star were obtained, showing a clear high- amplitude signal 
in the data. Though the number of measurements is small 
(and comparable with the number of free parameters to 
fit), a tentative Keplerian fit holds a period of 22.6 days, 
eccentricity of 0.50, and a semi-amplitude of 9.5 km s . 
The signal is compatible with a 92 Jupiter mass stellar 
companion to this 0.72 solar mass star. We note, however, 
that this orbital solution should be regarded with caution. 



4.2. Fast-rotating or giant stars 

HD221580: The 54 precise radial-velocity measurements of 
this star show a high-amplitude radial-velocity variation 
(total amplitude ^80ms _1 , rms=19.4ms~ 1 ). The data 
presents a rather complex behavior, and no clear period- 
icity could be found. Later detailed spectroscopic analy- 
sis has revealed that this star is a metal-poo r giant, with 
T eff =5322K, log g=2.68, and [Pe/H]=-1.13 (jSousa et all 
2010). The observed "noise" is thus likely to be related to 
intrinsic stellar variability. Very simila r stellar parameters 
were derived by iGratton et all (|2000i who classified this 
star as a red horizontal branc h giant. No reference to its 
evolutionary status is given in lNordstr5m et al.l (|2004). 



HD128575, HD197890, and BD- 00 4234: These three 
stars were excluded from the sample since the first (and 
only) obtained HARPS spectrum has shown that their 
CCFs were very large and shallow, or were very difficult 
to identify. This behavior is typical of fast-rotating stars. 
These objects are not suitable for a precise radial velocity 
search for planets. 



5. Results 

After removing the 16 targets discussed above, we are left 
with 88 stars that we consider suitable for a high-precision 
radial-velocity planet search. From these, however, 24 ob- 
jects have 5 or less measurements, making it difficult to 
make any firm conclusions. These are discussed separately 
from the others. For the remaining 64 stars (with 6 or more 
measurements), a more detailed analysis is done. 

5.1. Stars with <5 measurements 

Besides several of the cases discussed in Sect. @J 24 stars 
in our sample were observed only 5 or less times during 
the period of our survey (TablejSJ Fig.jS] one star, with 
one single measurement, was not included in the figure). 
With so few measurements, not much information can be 
extracted about these targets: a Keplerian fit, for instance, 
has 6 free parameters. We thus decided to separate all the 
targets with this reduced number of observations from the 
rest of the sample. 

In Table[5] we present the list of targets in such condi- 
tions, their weighted average radial velocity, the weighted 
rms of the measurements, the average error of the ob- 
tained radial velocities, the number of measurements, as 
well as their time span. A look at the dispersions and aver- 
age values shows that none of these stars seems to present 
strong radial-velocity variations. All dispersions are be- 
low 5ms -1 . Though the measurement time sampling may 
hide important signals, we can confidently exclude the ex- 
istence of giant planets in short-period orbits around these 
stars. A statistical approach to this issue will be taken in 
a separate paper. 

5.2. Tests for stars with at least 6 measurements 

For stars with 6 or more radial- velocity measurements (64 
objects - Figs. El El and|S]), we performed a series of sta- 
tistical tests, following the sa me methodology ap plied to 
the HARPS M-dwarf sample dBonfils et al.ll2010l First, 
to test whether the observed RVs (cr e ) are significantly in 
excess of the internal errors (<7j), we performed an F-test 
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Fig. 5. Radial velocity time series for stars with 2 to 5 radial velocity measurements. 
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Fig. 6. Radial velocity time series for stars with at least 6 radial velocity measurements. 



Santos et al.: The HARPS search for southern extrasolar planets 



1000 
JD (-2452500) 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 




1000 
JD (-2452500) 



600 800 1000 1200 1400 

JD (-2452500) 



1500 
_, 1000 

S" 500 

£ 

-500 
-1000 







I | I 1 I 1 | ! 1 , , | 1 I 1 J 3 


5 




10 


j i i i i 1 i i 

[\ 


1 1 | 1 1 1 1 | 1 . 








RV [m/s 
< a 


;1 1 . , - : 


1 
-5 


U 5 








, 1 , , , , 1 , , , , 1 , , , , : 


— D 




-10 


]~ T T T T T T 1 





500 1000 



1500 2000 
JD (-2452500) 



2500 3000 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 



1 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 



500 



1000 1500 
JD (-2452500) 



-10 - 1 

d 

500 



,.1 



1000 1500 
JD (-2452500) 



10 r . 

5 r J 
r 



1000 1500 2000 

JD (-2452500) 



1000 1500 2000 

JD (-2452500) 



$1 „ 




2500 500 



1000 1500 2000 

JD (-2452500) 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 



10 

— 5 - 

5 -5 b 
-10 



\ 1 ¥ 



1500 2000 
JD (-2452500) 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 



10 r 
— 5 r 
1 - 

-10 r " 



1000 1500 
JD (-2452500) 



1000 1500 
JD (-2452500) 



Fig. 7. Continuation of Fig.[BJ 
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Fig. 8. Continuation of Fig. [51 



and deriv ed the probability, P(F) to the F- value F= 



(see e.g. IZechmeister et al.ll2009i ). We also computed the 



Xconstant f° r the constant model as well as the probabil- 
ity of having xlonstant given the a,, P ' (xlonstant) ■ Th ese 
values are listed in Tables [6] and [7] For 16 out of the 64 
stars, both P(F) and P(Xconstant) are below 1%. These are 
the cases of CD-2310879, HD11397, HD59984, HD88725, 
HD107094,HD109310, HD123517, HD124785, HD144589, 
HD150177,HD171028,HD171587, HD181720, HD190984, 
HD215257, and HD224817. 

We note that the internal errors include both photon 
noise errors and errors from the wavelength calibration 
and instrumental drift correction. Furthermore, we added 
an instrumental error of 0.8 m s _1 to take the estimated in- 



trinsic precision of HARPS (guiding and internal errors) 
into account. Other sources of noise, like those coming 
from stellar intrinsic phenomena, are more difficult to es- 
timate, so were not considered here. 

Two tests were also done to estimate the probability 
that the observed radial velocities are best fitted by a lin- 
ear function of type rv = a x time + b. First, we used an 
F-test to derive the probability P(F s i ope ) that the use of an 
additional free parameter (2 instead of 1 in case of a con- 
stant fit) implies a statistically significant improvement 
when comparing X 2 s io P e to xlonstant- Second, we derived 
the false alarm probability using a bootstrap randomiza- 
tion approach. In this case, we generated random RV time 
series by shuffling (without repetition) the original radial 
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Name 


N 




cr e 


P(F) 


Xconstant 


P [~X.constant) 


Slope 


XsZope 


P ( Fa I ope ) 


FAP 






[ms- 1 ] 


[ms" 1 ] 








[m/s/yr] 








CD-2310879 


35 


4.2 


12.4 


< 10" y 


336 


< 10" y 


0.476 


336 


1.000 


0.961 


CD-571633 


7 


2.9 


2.9 


0.645 


5.3 


0.505 


-0.244 


4.5 


0.686 


0.107 


HD967 


34 


1.3 


1.8 


0.028 


56.9 


0.006 


-0.154 


55.4 


0.674 


0.009 


HD11397 


33 


1.4 


6.3 


< lO" 9 


462 


< lO" 9 


-4.579 


50.1 


< lO" 9 


< 1/720 


HD17548 


10 


2.8 


2.3 


0.823 


4.8 


0.851 


-0.244 


4.2 


0.496 


0.028 


HD17865 


21 


1.5 


1.6 


0.450 


27.4 


0.123 


-0.073 


27.3 


1.000 


0.548 


HD22879 


36 


1.3 


1.4 


0.468 


49.0 


0.059 


-0.297 


40.4 


< lO" 9 


0.077 


HD31128 


37 


3.7 


4.3 


0.245 


57.5 


0.013 


-0.547 


54.5 


0.035 


0.070 


HD40865 


30 


1.5 


2.3 


0.014 


48.3 


0.014 


0.062 


48.1 


1.000 


0.019 


HD51754 


21 


1.6 


1.9 


0.333 


29.8 


0.072 


-0.729 


22.4 


1.3 10 -4 


0.016 


HD56274 


14 


1.3 


2.3 


0.035 


35.8 


6.3 10 -4 


-0.783 


23.4 


0.002 


0.002 


HD59984 


45 


1.2 


2.4 


2.0 10 -6 


202 


< lO" 9 


0.712 


185 


6.0 10 -6 


0.019 


HD61902 


7 


2.9 


4.8 


0.217 


20.0 


0.003 


-1.362 


11.1 


0.106 


0.056 


HD68089 


7 


2.4 


2.6 


0.598 


7.2 


0.302 


-0.412 


6.6 


0.902 


0.586 


HD68284 


10 


1.9 


3.0 


0.167 


21.4 


0.011 


0.459 


17.6 


0.305 


0.026 


HD69611 


6 


1.8 


2.5 


0.392 


13.4 


0.020 


0.774 


10.4 


0.609 


0.365 


HD75745 


13 


2.8 


3.6 


0.252 


18.9 


0.092 


0.722 


12.9 


0.008 


0.023 


HD77110 


16 


1.4 


2.1 


0.089 


35.5 


0.002 


-0.964 


15.9 


2.0 10~ 6 


< 1/720 


HD78747 


26 


1.2 


1.6 


0.184 


40.5 


0.026 


-0.137 


39.3 


0.799 


0.099 


HD79601 


16 


1.6 


2.0 


0.282 


23.6 


0.071 


-0.569 


11.2 


4.0 10 -6 


< 1/720 


HD88474 


6 


3.9 


4.5 


0.534 


9.6 


0.088 


1.035 


5.9 


0.274 


0.177 


HD88725 


22 


1.2 


3.3 


1.3 10 -5 


179 


< lO" 9 


-1.879 


46.8 


< lO" 9 


< 1/720 


HD90422 


7 


2.6 


3.6 


0.347 


8.6 


0.199 


-0.318 


8.1 


0.970 


0.159 


HD91345 


8 


4.0 


4.7 


0.476 


13.1 


0.070 


-1.574 


4.8 


0.006 


0.021 


HD94444 


7 


2.3 


4.0 


0.181 


24.8 


3.7 10 -4 


-1.413 


13.7 


0.101 


0.162 


HD95860 


7 


2.3 


3.7 


0.204 


20.8 


0.002 


0.059 


20.3 


0.997 


0.900 


HD97320 


6 


3.1 


2.8 


0.726 


5.9 


0.318 


-0.283 


4.9 


0.741 


0.471 


HD97783 


6 


1.9 


2.5 


0.443 


15.3 


0.009 


-0.300 


13.2 


0.816 


0.748 



velocity data, preserving the same observing dates. On 
each obtained time series, we adjusted a slope and com- 
puted its x 2 value. The fraction of simulated data sets with 
X 2 lower than the observed one gives us the false-alarm 
probability (FAP) for the slope model. Thirteen (13) stars 
in our sample show both significant P(F s i ope ) and FAP val- 
ues: HD11397, HD56274, HD77110, HD79601, HD88725, 
HD107094,HD109310, HD113679, HD123517, HD144589, 
HD148211,HD190984, and HD215257. 

For stars with at least 12 measurements (37 stars), a 
test was also done to search for previously undetected pe- 
riodic signals in the data. For this we analyzed the general- 
ized lomb sc argle periodogram. We followed the p r escrip - 
tion used inlBonfils et~alT(l2010l). ICumming et all (|l999h . 
and Zechmeister fc Kiirsterl (|2009h . In brief, we created 
10 000 virtual time series by making permutations of the 
original radial velocity set. For each case we then com- 
puted a periodogram and located the highest power. We 
then derived the distribution of the power maxima. We 
finally attributed an FAP to the maximum power value 
found in the original data set by counting the fraction of 
the simulated power maxima that have a higher value. The 
resulting FAP values are listed in Table[5J together with 
the amplitude of the period and amplitude of the highest 
peak in the periodogram of the real data. For stars with 
known planetary companions, both FAP before and after 



the subtraction of the Keplerian solution are shown, while 
for stars presenting linear drifts the test was done only 
after the subtraction of the trend. This was done since a 
linear drift will introduce extra power at low frequencies, 
masking any shorter period significant peaks. Six (6) stars 
in our sample show the presence of significant peaks (FAP 
below or near 1%): CD-2310879, HD78747, HD107094, 
HD171028, HD181720, and HD190984. 

Finally, for all the targets a Keplerian analysis was 
done using yorbit (Segransan et al., in prep.), a code that 
uses an hybrid method based on a fast linear algorithm 
(Levenberg-Marquardt) and genetic operators (breeding, 
mutations, cross-over). This approach is important since 
in a periodogram an eccentric orbit is translated into not 
one but multiple peaks, making it difficult to identify the 
correct orbital parameters. 



5.3. Case-by-case analysis 

For 9 of the cases mentioned above, HD11397, HD123517, 
HD144589, and HD215257 (clear linear drifts), HD171028, 
HD181720, and HD190984 (previously announced plan- 
ets), CD-2310879 (active star), and HD107094 (candi- 
date new planet) , the cause of the observed variation will 
be discussed in Sects. l5XTl 15X21 and |5"X1 A dis- 
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Table 7. Variability tests for stars with at least 6 measurements (HD numbers above 100 000). 



Name 


N 


Oi 




P(F) 


2 

Xconstant 


P * {.Xconstant ) 


Slope 


'2 

Aslope 




FAP 






[ms x ] 


r —It 

[ms x \ 








[m/s/yr] 








HD102200 


7 


5.4 


8.1 


0.279 


7.7 


0.260 


0.835 


4.5 


0.127 


0.023 


HD 104800 


6 


1.8 


1.5 


0.795 


3.2 


0.663 


-0.677 


2.1 


0.325 


0.129 


HD107094 


14 


1.8 


599 


< 10 -9 


1.8 10+ 6 


< 10 -9 


429.034 


1.3 10+ 4 


< 10 -9 


< 1/720 


HD108564 


6 


1.6 


1.3 


0.791 


3.8 


0.573 


-0.436 


2.5 


0.321 


0.176 


HD109310 


15 


1.6 


3.7 


0.003 


78.5 


< 10~ 9 


-1.460 


47.9 


3.1 10~ 4 


0.003 


HD109684 


6 


1.9 


2.5 


0.450 


8.0 


0.158 


-0.202 


6.3 


0.645 


0.170 


HD111777 


6 


1.5 


1.2 


0.814 


4.2 


0.522 


-0.355 


2.9 


0.403 


0.299 


HD113679 


6 


2.2 


2.7 


0.496 


9.8 


0.082 


1.697 


0.4 


1.7 10~ 4 


< 1/720 


HD123517 


9 


3.1 


49.5 


< 10 -9 


2039 


< 10" 9 


34.103 


3.2 


< 10 -9 


< 1/720 


HD124785 


17 


2.2 


4.7 


0.004 


78.8 


< 10" 9 


-1.794 


68.3 


0.080 


0.100 


HD126681 


13 


2.4 


3.7 


0.116 


23.6 


0.023 


-0.788 


20.1 


0.201 


0.057 


HD126793 


7 


2.3 


4.1 


0.169 


14.1 


0.029 


0.737 


7.7 


0.100 


0.020 


HD126803 


7 


2.6 


4.8 


0.143 


12.9 


0.044 


0.414 


8.7 


0.252 


0.022 


HD 134440 


10 


2.4 


2.5 


0.562 


9.7 


0.375 


-0.171 


9.0 


0.857 


0.083 


HD 144589 


11 


4.0 


34.7 


< 10 


604 


< 10 


25.404 


14.2 


< 10 


0.001 


HD148211 


31 


1.6 


2.5 


0.011 


76.0 


7.0 10~ 6 


-0.600 


67.8 


8.2 10~ 4 


0.006 


HD148816 


6 


1.6 


1.8 


0.569 


5.7 


0.341 


0.198 


3.6 


0.305 


0.070 


HD 149747 


9 


4.7 


4.6 


0.666 


10.4 


0.237 


-0.904 


7.9 


0.213 


0.142 


HD150177 


30 


1.6 


4.8 


< 10~ 9 


327 


< 10~ 9 


-0.091 


326 


1.000 


0.313 


HD167300 


9 


1.8 


2.1 


0.479 


12.4 


0.132 


-0.135 


11.9 


0.968 


0.596 


HD171028 


48 


1.5 


33.8 


< 10~ 9 


2.9 10+ 4 


< 10~ 9 


8.624 


2.7 10 +4 


1.0 10~ 6 


0.208 


HD171587 


14 


1.3 


6.0 


2.0 10~ 6 


261 


< 10 -9 


1.791 


225 


0.185 


0.175 


HD175607 


7 


1.7 


1.2 


0.897 


4.4 


0.618 


1.082 


1.5 


0.015 


0.051 


HD181720 


29 


1.8 


6.3 


< 10 -9 


617 


< 10 -9 


-0.899 


587 


0.247 


0.107 


HD190984 


46 


2.3 


39.1 


< 10" 9 


2.0 10+ 4 


< 10" 9 


20.284 


2910 


< 10" 9 


< 1/720 


HD197083 


12 


1.6 


2.7 


0.080 


45.1 


5.0 10~ 8 


1.517 


30.8 


0.016 


0.278 


HD197197 


21 


1.7 


2.9 


0.022 


49.0 


3.1 10 -4 


-0.609 


46.9 


0.721 


0.061 


HD199288 


14 


1.0 


1.0 


0.509 


15.7 


0.264 


-0.373 


11.4 


0.009 


0.064 


HD199604 


6 


1.8 


1.8 


0.707 


3.6 


0.610 


2.391 


1.2 


0.041 


0.055 


HD199847 


7 


2.0 


2.9 


0.331 


9.3 


0.156 


1.331 


4.5 


0.059 


0.021 


HD206998 


6 


2.5 


2.8 


0.563 


5.7 


0.332 


-1.071 


2.0 


0.047 


0.018 


HD207869 


17 


1.4 


1.4 


0.613 


19.1 


0.266 


-0.305 


17.4 


0.339 


0.373 


HD210752 


17 


1.4 


2.1 


0.104 


39.0 


0.001 


0.378 


34.3 


0.126 


0.039 


HD215257 


37 


1.6 


14.4 


< lO" 9 


1870 


< lO" 9 


-15.045 


343 


< lO" 9 


< 1/720 


HD218504 


15 


1.3 


2.0 


0.122 


30.6 


0.006 


-0.449 


26.4 


0.140 


0.095 


HD224817 


30 


1.4 


3.4 


6.0 10~ 6 


109 


< lO" 9 


-0.134 


105 


0.691 


0.090 



cussion of the remaining stars presenting some sort of vari- 
ability according to the statistical tests is done below. 

HD56274 (Nmes=14): The radial velocities of this 
star suggest the existence of a very low-amplitude drift 
(~— 0.8m/s/yr), and our data does not allow us to reach 
any firm conclusions. No other significant signal is found 
in the data. No planet has been detected around this 
star by other high-prec ision radial velocity surveys (e.g. 
Fischer fc ValentfeOOFilh 

HD59984 (Nmes=45): This star shows an rms that 
is clearly above the average photon noise error. A pe- 
riod search (using both a periodogram and a Keplerian 
analysis) did not find any significant periodic signal in 
the data. The high dispersion observed may likely be ex- 
plained by t he ey olutionary status of the star. Indeed, 
Sousa et al. I (l2010h derive a surface gravity of 4.18 dex for 
this T e// =5962K metal-poor star ([Fe/H]=-0.69), indi- 
cating that it may be slightly evolved out of the main 
sequence. Evolved stars are known to have higher os- 



cillation and granulation noise (e.g. Mayor et al 



Christensen-Dalsgaardll2004l Dumusaue et al. 201C). 



2003 



HD77110 (Nmes=16): This star shows a tentative 
very low-amplitude drift in radial velocity (amplitude 
~— 1 m/s/yr). Such trend is not discussed in the litera- 
ture, although to our knowledge this star has not been in- 
cluded in any other high-precision radial- velocity survey. 
No significant periodic signal is found in the data. 

HD78747 (Nmes=26): The analysis of the generalized 
Lomb Scargle periodogram (TablejU) suggests the presence 
of a 3.1 day period signal in the data. It is indeed possible 
to find a satisfactory Keplerian fit with this period, an 
amplitude of 2.0ms _1 , and eccentricity of 0.26. In order 
to confirm this signal we obtained a series of 14 new radial 
velocities from April to July 2010. The analysis of the data 
(to be published in a separate paper) does not reveal any 
3-day signal. 

HD79601 (Nmes=16): This is another case of low- 
amplitude drift (—0.6 m/s/yr), whose cause is not identi- 
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Table 5. Stars with 5 or less measurements. 



Table 8. Periodogram variability tests for stars with at 
least 12 measurements. 



Star 


<RV>t o-fRVH 


^cr mes ^> 


N 


Span 
[days] 










kms" 1 ] 


ms" 1 ] 


[ms- 1 ] 




Star 1 " 


FAP 


Period 


Amplitude 


BD-084501 


84.469 


1.7 


7.4 


3 


1077 




[%] 


[days] 


[ms- 1 ] 


BD+062932 - 


144.809 


1.0 


2.2 


4 


1036 


CD-2310879 


1.08 


4.63 


14.2 


BD+063077 


46.843 


- 


7.3 


1 


— 


HD967 


41.09 


2.48 


1.3 


BD+083095 


-51.148 


3.2 


3.3 


3 


1378 


HD11397 d 


91.19 


2.50 


1.2 


CD-451246 


108.847 


3.8 


6.3 


4 


1001 


HD 17865 


22.16 


4.30 


1.8 


HD38510 


183.828 


2.6 


2.8 


5 


1502 


HD22879 


7.53 


6.11 


1.4 


HD105004 


121.495 


3.9 


3.1 


5 


1084 


HD31128 


29.77 


16.81 


4.0 


HD111515 


2.657 


0.9 


2.0 


5 


1458 


HD40865 


76.46 


8.69 


1.3 


HD 119949 


-39.327 


2.6 


3.0 


5 


1463 


HD51754 


91.65 


18.96 


1.6 


HD121004 


245.318 


1.7 


2.0 


5 


1470 


HD56274 d 


90.40 


4.78 


1.9 


HD128340 


2.997 


2.0 


2.7 


5 


1031 


HD59984 


3.80 


10.08 


2.5 


HD129229 


-15.932 


1.9 


2.0 


5 


1031 


HD75745 


28.95 


307.83 


3.0 


HD131653 


-67.090 


1.7 


1.2 


4 


1375 


HD77110 d 


92.28 


7.06 


1.4 


HD 134088 


-61.111 


1.0 


1.3 


4 


1377 


HD78747 


0.94 


3.09 


1.7 


HD145344 


-11.875 


1.8 


3.6 


5 


1378 


HD79601 d 


38.62 


9.56 


1.1 


HD145417 


8.799 


1.3 


2.2 


5 


1378 


HD88725 d 


54.37 


3.44 


1.5 


HD147518 


-42.127 


0.6 


2.0 


4 


1373 


HD107094 d 


<io- 4 


1898.28 


80.3 


HD193901 


171.416 


2.2 


2.8 


3 


1371 


HD109310 d 


5.93 


77.34 


6.0 


HD195633 


-45.832 


3.1 


3.0 


4 


346 


HD 124785 


55.12 


41.84 


6.6 


HD196892 


-34.498 


4.1 


2.5 


3 


602 


HD126681 


17.71 


1034.04 


6.2 


HD197536 


-17.117 


1.4 


2.0 


3 


1286 


HD148211 


53.78 


63.39 


1.9 


HD199289 


-6.144 


1.2 


2.6 


5 


1123 


HD 150 177 


2.46 


147.62 


5.6 


HD207190 


-11.314 


1.3 


2.5 


5 


1365 


HD171028 


<10" 4 


500.55 


47.1 


HD223854 


19.596 


2.4 


1.8 


4 


1401 


HD171028 P 


3.81 


16.71 


2.3 


f Weighted values 












HD171587 


71.19 


4.30 


5.1 












HD181720 
HD181720 P 
HD 190984 


<10~ 4 
50.68 
<10 4 


973.47 
56.25 
1948.73 


8.9 
1.3 
61.6 


fied. No significant periodic 


signal was found 


in 


the data, 


HD190984 P 


36.50 


33.29 


3.1 


even after subtracting the linear drift. No reference for a 


HD 197083 


21.06 


87.60 


7.1 


possible companion was found in 


the literature. 




HD197197 


9.67 


5.36 


2.8 


HD88725 (N 


mes=22): The r 


adial velocity 


measure- 


HD199288 
HD207869 


77.67 
26.90 


59.37 
15.67 


2.0 
2.5 


ments of this star present 


a low 


amplitude linear trend 


HD210752 


25.38 


2.60 


1.9 


(— 1.9m/s/yr). The analysis of the residuals of a linear fit 


HD215257 d 


20.28 


8.05 


3.7 


to the data show the presence of a peak around 3 days. 


HD218504 


80.32 


4.96 


1.9 


However, the FAP of this peak 


suggests that 


it is not 


HD224817 


8.62 


9.15 


2.0 



significant. The small number of measurements precludes 
any further conclusions. No planet was found around this 
star in other high pre cision planet search surveys (e.g. 
Fischer fc Valentill2005l) . 

HD109310 (Nmes=15): A low-amplitude drift is sug- 
gested from our radial velocity data (— 1.5 m/s/yr). The 
residuals of a linear fit to the data show an rms of 



, still above the average radial velocity error 
. A frequency analysis of the residuals does 



~2.5ms" 
of 1.6ms 

not reveal, however, any significant signal. Having so few 
measurements does not allow us to make any further con- 
clusions about this case. 

HD113679 (Nmes=6): A long-term and low-amplitude 
(~1.7ms _1 ) linear trend is observed for this star over a 
span of ^1500 days. The reduced number of data points 
obtained does not allow reaching any firm conclusions 
about the observed variation. We found no reference in th e 
literature to any variability (e.g. iNordstrom et al. 20041 ). 



although this star was not included (to our knowledge) in 
any other high-precision radial-velocity survey. 



* Names with suffix d and p denote the results after subtraction 
of the linear drift and Keplerian solutions, respectively. 



HD124785 (Nmes=17): This star presents a clear 
excess rms in the radial velocities, when compared to 
the average error. As for HD59984, however, the atmo- 
spheric parameters derived by Sousa et al. (T e y/=5867K, 
log g=4.20) suggest that this star has already evolved off 
the main sequence. This may explain the excess power 
observed, and is supported by not finding any significant 
periodicity in the data. 

HD148211 (Nmes=31): The statistical tests done 
above suggest there is a significant drift with a very 
small amplitude (-0.6ms -1 ). This signal is visually very 
marginal, though. No significant periodicity is found in 
the data. Our statistical tests also do not indicate any 
significant radial- velocity excess. 

HD150177 (Nmes=30): A high level of residuals is seen 
in the data for this star. No significant periodic signal(s) 
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Fig. 9. Radial velocity time series for the 3 stars with 
planets discussed in Sect. 15.3. T1 



was found though. The early spectral type and evolution- 
ary status suggested by the precise atmospheric param- 



eters derived by Sousa et al. (T e j/=6216K, logp=4.18) 
indicate that the observed noise may be due to stellar os- 
cillation and granulation phenomena. 

HD171587 (Nmes=14): The radial velocity measure- 
ments of this star show a clear radial-velocity rms ex- 
cess with respect to the measurement errors. The analysis 
of the data show that, excluding the point at JD=54638 
(~20ms _1 above the average of the other measurements), 
a clear periodic signal appears that is best fit with a 
Kcplerian function with P=4.49days, K=4.0ms _1 , and 
e=0.57. The residuals of this fit are only 0.56ms _1 (com- 
pared to the average error of 1.3 ms -1 ), making us sus- 
pect that there are not enough points to obtain a reli- 
able fit. This signal would be compatible with the pres- 
ence of a 7.1 Earth mass planet orbiting this 0.76 Mq star 
(T e// =5412K, logg=4.59, [Fe/H]=-0.64). The analysis 
of the Ca n H and K lines suggests that this star is chromo- 
spherically very active (log R' HK =— 4.75). Interestingly, 
this activity value indicates a rotational period of 21 days 
( Noves et al.lll984 ). ~5 times longer than the observed sig- 
nal. A periodogram of the individual Smw values shows 
a peak around 23.5 days, with a significance level of 93%. 
This peak may be the signature of the rotational period 
of the star. We note, however, that the Keplerian fit de- 
scribed is only based on 13 points (excluding the above- 
mentioned date), and the time sampling is very wide (in- 
dividual points are very separated in time) . A series of ad- 
ditional data points were obtained in 5 consecutive nights 
in July 2010 under ESO program 085.C-0063 (to be pub- 
lished in a separate paper). The results do not show any 
clear 4.5 day periodicity. 

HD224817 (Nmes=30): According to our statisti- 
cal tests, this star shows an rms above the expected 
value when taking the average internal errors into ac- 
count. However, this is largely due to the only point at 
JD=53289, without which no significant variation is seen. 

5.3.1. Previously announced planets 

Three stars in our sample were previ ously announced as 



2007 



harb oring long-period giant planets (pantos et al 
2010h . see Fig.EJ For one of these, new HARPS radial 



velocities were obtained after the discovery paper. 

HD 171 028: This star was the first in our sample to be 



anno unced to have an orbiting giant planet (pantos et al 



20071 ). Since the announcement paper, 29 new radial veloc- 
ity measurements were obtained (in a total of 48), allowing 
us to put stronger constraints to the orbital parameters of 
HD171028b: period of 550 days, eccentricity of 0.59, and 
minimum mass of 1.98 (supposing a stellar mass of 



iij.ij.iiiiinj.il iiicjoo uj j-.^yu ,J wp V l -i}-'^ wc 'iii t 3 ci □ t^uaj iiicioo UJ 

1.01 M - ISousa et aHl2010l) . In TableEXUwe present our 



solution. A Lomb Scargle Periodogram of the residuals re- 
veal the presence of a non significant peak at ^17 day 
period. This may be related to the rotational period of 
the star. We tried to fit a Keplerian function to the resid- 
uals of the Keplerian fit, but no satisfactory solution was 
found. 
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Table 9. Elements of the fitted orbits for known planets. 
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HD 171028 


HD 181720+ 


HD190984 




p 


550±3 


956±14 


4885±1600 


[d] 


T 


2 454187±1 


2 453 631±30 


2 449 572±1600 


[d] 


a 


1.32 


1.78 


5.5 


[AU] 


e 


0.59±0.01 


0.26±0.06 


0.57±0.10 




V r 


13.641±0.001 


-45.3352±0.0004 


20.269±0.004 


[kms -1 ] 


OJ 


304±1 


177±12 


318±5 


[degr] 


K\ 


60.6±1.0 


8.4±0.4 


48±1 


[ms- 1 ] 


him) 


6.61 1(T 9 


0.053 10 -9 


30.95 


[M s ] 


a{0 - C) 


2.3 


1.37 


3.44 


[ms" 1 ] 


N 


48 


28 


47 




m2 sini 


1.98 


0.37 


3.1 


[M Jup ] 



f Orbital parameters from Santos et al '. (2010) 



HD 181720 and HD 190984: giant planets were found 
around thes e two stars and recently announced in 
Santos et al. (|2010h . The orbital parameters for the two 
planets are listed in TableESHJ We refer to the discovery 
paper for more details. 

5.3.2. High-amplitude linear drifts 

A few stars in our sample present very clear long-term 
trends with total amplitudes of at least 10 to 20 ms -1 . 
These are compatible with the existence of long period 
giant planets or brown dwarfs. However, the insufficient 
time coverage of the data does not allow making any con- 
clusion about the origin of the observed signal. 

HD11397: The 33 radial velocities of this star show the 
signature of a clear long-term trend with a total amplitude 
of 33 ms -1 over 2075 days. When excluding the first two 
measurements, the trend shows a slight curvature. The 
data could be well fit by a linear slope, together with 
a low-amplitude long-period signal. However, the lower 
quality of the first two data points (photon noise error 
of 2-4 ms -1 and obtained with short exposure times not 
allowing supression of stellar oscillation noise) keeps us 
from deeper insight into the observed signal. No signifi- 
cant short-period signal is found in the data. A quadratic 
fit that excludes the first two points shows a low rms of 
1.27ms -1 , not strongly above the average photon noise 
(0.90ms -1 ). More data are needed to clarify this case. 

HD123517: The 9 radial velocity measurements of this 
star reveal a clear long-term linear drift (~150ms -1 over 
~1600days). The residuals of a linear fit to the data 
(1.5ms -1 ) are smaller than the average photon noise error 
of the measurements (2.3ms -1 ; this is one of the faintest 
stars in our sample with mv=9.6), suggesting that no de- 
tectable short period signal is hidden in the data. More 
data are needed to confirm the source of the observed ra- 
dial velocity variation. 

HD144589: We obtained 11 radial velocity measure- 
ments of this star, spanning a total of ^1500 days. The 
data shows the presence of a clear linear trend, with a fit- 
ted slope of 25 ms -1 yr -1 . The residuals of a linear fit to 



the data (4.0 ms- 1 ) are slightly above the average photon 
noise of the measurements (3.4ms- 1 ). The stellar atmo- 
spheric parameters derived f or this star ( T e ff= 6372K, 
log£f=4.28, [Fe/H]=-0.05) bv lSousa et aH (<2010l ) suggest 
that this is one of the hottest stars in our sample and likely 
slightly evolved off the main sequence. Stellar oscillations, 
along with the relatively high projected rotational velocity 
(5.6 kms -1 - derived from the HARPS cross-correlation 
function), may explain the excess observed scatter. 

HD215257: The 37 radial velocities of this star (av- 
erage error of 1.6 ms- 1 ) indicate the presence of a low 
amplitude linear drift (— 15ms -1 yr -1 over 1700 days). A 



signal with a similar slope was also found bv lSozzetti et al 



(2009) using lower precision data from their Keck-HIRES 
survey. The high residuals (3.98 ms- 1 ) around the linear 
trend indicate an extra signal in the data. No reliable 
Keplerian fit was found, though, and no significant peak is 
found in a generalized Lomb Scargle periodogram. Since 
the star is an early-G/late-F dwarf with T e ff =6052 K, 
log 3=4.46, and [Fe/H]=-0.63 dSousa et al.ll2010h . at least 
part of this signal may come from stellar oscillation and 
g ranulation "noise" ( more important i n ear ly type stars 
Mayor et all l2003t iDumusaue et alj|2010l) . More data 



are needed to confirm that this is the only source of the 
observed rms. 



5.3.3. Activity-induced periodic signals 

CD- 23 10879: The 35 radial- velocity measurements of 
this F dwarf (T e// =6788K, log 5 =4.67, [Fe/H]=-0.24) 
present an rms of 12.9 ms -1 , clearly above the median 
photon noise (3.9ms -1 ). A Lomb Scargle periodogram 
shows a significant peak around 4.5 days, well fit by a 
Keplerian function (see Fig. [T0|) with P=4. 626 days, eccen- 
tricity of 0.092, and semi-amplitude of 14.4 ms- 1 (resid- 
uals of 3.96ms- 1 , similar to the photon noise). This sig- 
nal is compatible with the existence of a planet with 42 
times the mass of the Earth orbiting this 1.20 Mq star. 
However, a Lomb Scargle periodogram of the bisector in- 
verse slope (BIS) of the HARPS cross-correlation func- 
tion shows a similar peak at ^4.5 days. A clear negative 
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Fig. 10. Top: Phase folded radial-velocity measurements 
of CD-2310879 with a period of 4.626 days. The line rep- 
resents the best Keplerian function fit. Bottom: Radial 
velocities as a function of the HARPS cross-correlation 
bisector inverse slope. The line represents a linear fit to 
the data. 

correlation is also seen between BIS and the radial veloc- 
ities (slope=-0.94±0.23 - Fig. [TO]), typical of the signal 
expected from activity-in duced radial-velocity variations 
(see case of HD 166435 



Queloz et al.l l2000h . The high 



value of the activity index obtained from the HARPS spec- 
tra, log R' HK —— 4.60, suggests a rotational period around 
4 days, similar to t he observed radial-velocity periodicity 
(jNoves et al.lll984f ). We thus conclude that the best ex- 
planation for the observed radial-velocity variation of this 
star is given by active regions modulated by stellar rota- 
tion. 

5.3.4. Ambiguous cases: new candidate planets 

HD 107094'- A total of 14 radial velocities were obtained for 
this star, showing a clear long-term (quasi-linear) radial- 
velocity signal. Twelve of these measurements were ob- 
tained within the HARPS GTO program, and a last two 
in July 2010 under ESO program 085.C-0063. In a first 
analysis, the total amplitude (~2.5kms _1 over 2355 days) 
suggests there is a stellar or sub-stellar companion or- 
biting this 0.8 solar mass star. A single Keplerian fit 
to the data shows that it is well fit by a long-period 
Keplerian function, though the orbital parameters are dif- 
ficult to constrain (see Fig.QT]). The residuals of this fit, 
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Fig. 11. Top: Radial velocity measurements of HD107094 
and best 1-Keplerian fit to the data. Middle: same as above 
but after fitting a linear trend, together with a Keplerian 
function. In both cases the residuals are shown. Bottom: 
phase-folded radial velocity measurements with the period 
of the candidate planet. 



on the order of 2.1ms _1 , are still high, and they called 
attention to a different possible orbital solution to the 
system. Indeed, a fit of the linear drift together with 
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Table 10. Elements of the fitted orbit for HD107094. 



p 


1870±34 


[d] 

L 1 


T 


2 454 765±29 


[d] 


a 


2.76 


[AU] 


e 


0.13±0.02 




V r 


15.33±0.05 


[kms" 1 ] 


CJ 


46±3 


[degr] 
[ms" 1 ] 


Ki 


88±5 


Mm) 


128.3 10" 9 


[M ] 


o(0 - C) 


1.10 


[ms' 1 ] 


N 


14 




iri2 sini 


4.5 


[M Jup ] 



a Keplerian function reveals the presence of a possible 
shorter period solution, with P^1870days, K=88ms _1 , 
and e=0.13 (Table lO^ . This signal, first detected us- 
ing only the GTO data, is compatible with the presence 
of a 4.5 Jupiter mass companion. The residuals of this 
fit (1.10 ms -1 ) are significantly better than in the previ- 
ous case. The fitted linear trend (which, to our knowl- 
edge, has not been previously reported in the litera- 
ture) has a slope of 402±lm/s/yr. More data are now 
needed to disentangle the two possible solutions, but we 
are likely seein g a new giant pla net. Interestingly, with 
[Fe/H]=-0.51 (ISousa et alJl2010h . HD107094 is close to 
the high-metallicity limit of our sample. 

HD197083: Although the statistical analysis done for 
this star does not indicate any strong signal, a visual in- 
spection of the radial velocities hints at the presence of a 
long-term periodic signal (P~1500-2000days), see Fig. [5] 
The number of points covering the first ^1000 days is 
very reduced (3) with high error bars. More data are thus 
clearly needed to confirm this signal. 

6. Discussion and conclusions 

In this paper we present the overall results of our HARPS 
program to search for planets orbiting a sample of metal- 
poor stars. In a total of 104 stars, 3 giant planets were 
confirmed, and one more very promising giant planet can- 
didate is presented. These discoveries significantly increase 
the number of known giant planet s orbiting metal-poo r 
stars (see also case of HD 155358 - ICochran et~"aTll2007l) . 
Several binary stars were also found in the sample. 

A proper statistical analysis of our results requires not 
only identifying the planetary signals in our time series but 
also assessing their detection limit; i.e. what planets, as a 
function of mass and period, can be rejected considering 
our observations. Although this is the topic of a forthcom- 
ing paper, partial figures can already be given. Excluding 
the 16 stars discussed in Sect. SI (targets not suitable for a 
high-precision radial- velocity planet search), our sample is 
constituted of 64/87 stars with at least 6/3 radial velocity 
measurements. We have detected at least 3 giant planets 
with masses greater than 0.3 Jupiter masses and with pe- 
riods ranging from a few hundred to a few thousand days. 
Of course, time series with only a few measurements (~3- 



6) cannot always rule out a Jupiter-mass companion on a 
short or moderate period orbit. Moreover, we have a pos- 
sible fourth detection and several more stars with a lin- 
ear RV drift indicative of a long-period companion. Three 
(3) detections over 87 targets, or SAtHm is therefore 
a lower estimate to the frequency of Jupiter-mass planets 
orbiting metal-poor stars. If we only consider stars with 
[Fe/H] between —0.40 and — 0.60dex (an interval where 
the stars with planets are in our sample), then only 34 
stars with at least 3 measurements exist, and the percent- 
ages increase to 11. 3^3%. 

We can additionally note that no planet has been de- 
tected around the 32 stars with [Fe/H]< — 0.60dex and 
at least 6 measurements. For these stars it is reasonable 
to assume completeness for Jupiter-mass planets on short 
or moderate period orbits (< 100 days). Those stars have 
metallicities from —1.4 to — 0.6dex, with a median value 
of — 0.69dex. A null detection for that subsample implies 
a maximum occurrence of Jupiter-mass planets (on short 
or moderate orb i tal pe r iods) o f ~5%. 



Santos et aD (|200lL l2004tJl have shown that the fre- 



quency of giant planets orbiting solar type dwarfs is a 
strong function of the stellar metallicit y. This trend, con- 
firmed by several other studies (e.g. iFischer fc Valenti 
120051 I Johnson et aHl2010h . makes these facts surprising. 
The above studies have also shown that the frequency of 
giant planets orbiting star s with metallicity near — 0.5dex 
is ^3% (see also review by lUdrv fc Santosll2007l ). strongly 
below the ~11% mentioned above (though both results 
are still likely compatible within the error bars) . We have 
to take into account, however, that the exquisite preci- 
sion of HARPS has increased the planet detection rate 
when compared to the results of the lower precision pro- 
grams on which the existence of the metallicity-giant 



planet correlation has been based (pantos et alj 12004b 



Fischer fc Valentil 120051 ). The present results, although 
based on small number statistics, may hint at a much 
higher frequency of giant planets around solar type stars 
(of all metallicities). 

The stars with giant planets discovered in our 
survey (HD 171028, HD 181720, HD 190984, and likely 
HD 107094), all have metallicity values on the high- 
metallicity side of our sample (—0.48, —0.53, —0.49, and 
—0.51, respectively), see Fig.Q] This suggests that even 
for metal-poor stars the frequency of giant planets is a 
rising function of the stellar metallicity, a result that does 
not confirm former suspicions that the metallicity-giant 
planet correlation could b e flat for [Fe/H] values below 
solar ( Santos et al. 2004b). This result further suggests 



that the core accretion model of planet formation is still 
at play in these low m etallicity values (jlda fc Lin||2004b ; 
Mordasini et al.ll2009aT) . Since no planets were discovered 
in our sample at [Fe/H] significantly below — 0.5dex, and 
only one giant planet wa s found around a star s with lower 
metallicity (— 0.68dex - ICochran et al.ll2007t ). this could 



4 The 1-sigma error bars were obtained using a binomial dis- 
tribution. 
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also hint that we are close to the metallicity limit below 
which no giant planets can be formed. More data are, how- 
ever, needed to settle this issue. 

As mentioned in Sect. [21 it is interesting to ob- 
serve that two of the stars discussed in Sect. 15.3.21 the 
ones that present long-term radial-velocity linear drifts 
(HD 123517 and HD 144589), have the highest metallic- 
ity values (+0.09 and —0.05 dex) among all the objects in 
our sample (see also Fig.[lJ. Both present relatively low- 
amplitude linear trends, suggestive of the presence of long- 
period giant planets. If such detections are confirmed, this 
result would be in line with the metallicity-giant planet 
correlation. 

Interestingly, no short period planet was found in our 
data. All the 3 (or 4) giant planets discovered orbit in 
long-period orbits. This suggests that giant planets in 
short-period orbits are not common around stars with low 



metal licity, in agreement with the results of lSozzetti et al 



(2009). The lack of such detections is, however, expected, 
since the frequency of hot .jupiters in rad ial velocity sur- 
veys is only ~1% (|Udrv fc Santos! l2007h . Only 87 stars 
were effectively monitored in the present survey, preclud- 
ing a significant statistical meaning for this lack of detec- 
tion. 

Theoretical studies (e.g. iMordasini et al 



2009a), 



backed up by recent observations (e.g. ISousa et al J 2008), 
suggest that the frequency of low-mass planets may in- 
crease for stars in the metallicity range of our sample when 
compared to solar metallicity objects. That no Neptune or 
super-Earth mass planet has been discovered up to now 
in our sample does not contradict this conclusion. Instead, 
the observing strategy used (data sampling and exposure 
times) was simply not adequate for the purpose of find- 
ing such low-mass objects. A program with HARPS is 
presently being done to fill this gap. 
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